Colloidal motility and pattern formation under rectified diffusiophoresis 
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In this letter, we characterize experimentally the diffusiophoretic motion of colloids and A- 
DNA toward higher concentration of solutes, using microfluidic technology to build spatially- and 
temporally-controlled concentration gradients. We then demonstrate that segregation and spatial 
patterning of the particles can be achieved from temporal variations of the solute concentration 
profile. This segregation takes the form of a strong trapping potential, stemming from an osmot- 
ically induced rectification mechanism of the solute time-dependent variations. Depending on the 
spatial and temporal symmetry of the solute signal, localization patterns with various shapes can be 
achieved. These results highlight the role of solute contrasts in out-of-equilibrium processes occuring 
in soft matter. 

PACS numbers: 



Diffusiophoresis is the mechanism by which particles 
and molecules drift along solute concentration gradients 
[TH2]. It belongs to the general class of phoretic trans- 
port phenomena, such as electro- and thermo- phoresis 
[4]-[8], i.e. the migration of particles under gradients of 
an external thermodynamic variable (electric potential, 
temperature, ...). Physically, the diffusio-phoretic migra- 
tion under solute gradients takes its origin in an osmotic 
pressure gradient occurring within the diffuse interface at 
the surface of the particle [TJ [2j |9] : this leads to a drift ve- 
locity proportional to the solute gradient Vc, in the form 
Vdp = /iVc, with /i a mobility. The diffusiophoretic 
transport was recognized since the pioneering work of 
Derjaguin and later by Anderson and Prieve [1 . How- 
ever its implication in out-of-equilibrium phenomena oc- 
curing in soft matter systems - in which solute contrasts 
are ubiquituous - has been barely explored up to now. 
Several works recently revealed the role played by trans- 
port induced by solute contrasts in out-of-equilibrium 
processes. It was e.g. shown to be at the origin of the 
strongly enhanced effective diffusion of colloids under so- 
lute gradients [2 with potential implications in microflu- 
idics [10 . Its possible interplay with other transport phe- 
nomena was also demonstrated, concuring for example to 
induce a sign reversal of thermophoretic transport of col- 
loids in the presence of polymers [6 j. Its role in skin 
formation during evaporation of mixtures was suggested 
[TT] . Finally, in the context of the efforts to develop 
autonomous microsystems, diffusiophoresis was used to 
harness chemical power to produce self-propelled artifi- 
cial swimmers [12], using chemical gradients as driving 
forces. 

In this paper we pursue the exploration of this phe- 
nomenon and demonstrate segregation and pattern for- 
mation of colloids and macromolecules on the basis of 
diffusiophoretic transport. We evidence a localization 
process of particles -here, colloids and A-DNA-, tak- 
ing its origin in the rectification of the underlying time- 
dependent concentration variations of a solute specie. 
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FIG. 1: (a) Experimental setup (b)-(Left) A 3-channel and 
gel matrix setup, superimposed on fluorescence intensity im- 
age measured for the stationary gradient of fluorescein with 
Ci-C 2 = 10~ 4 M O 300/xm, t 800 fim, scale bar 300/xm). 
(b)- (Right) Fluorescence intensity profiles measured in sta- 
tionary state, showing the expected linear profile. Actuation 
of the microfluidic switch allows to invert the concentration 
gradient. 



The shape of the pattern is shown to be strongly de- 
pendent on the temporal symmetry of the underlying 
solute concentration signal. A theoretical model based 
on a Smoluchowski description allows to reproduce the 
experimental observations. 

The thorough exploration of diffusiophoresis requires 
the building-up of controlled salt concentration gradients. 
To this end, we have developed a microfluidic experimen- 
tal setup sketched in Figjl] inspired from Ref. IT3] . This 
set-up allows to impose stationary, as well as temporally 
swit enable, solute gradients. A three channels device is 
molded in agarose gel, allowing free diffusion of solute 
without convective flow. A double syringe pump with 
two solutions of salt of concentration C\ and C2 fills the 
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FIG. 2: Diffusiophoretic transport of fluorescent colloids and 
A-DNA under a LiCl gradient (AC s [LiCl] = \C 2 - Ci| = 
lOOmM. £ = 800/rni, scale bar 100/xm). (a)-(c)). Motion of 
particles under a salt gradient, sketched by the lateral bar, 
towards higher salt concentrations. Images separated by 90s 
for the colloids and images at t =100, 150, 200, 300s for A- 
DNA. (b)-(d) Time evolution of the particles population lo- 
cation. Experimental data (symbols) are fitted according to 
the theoretical description (dashed lines), with Ddp as the 
only fitting parameter (see text). In (b), open symbols corre- 
spond to a subsequent migration, and fully superimposes on 
the previous results. The solid straight line in (b) is a guide 
line corresponding to constant drift velocity. Shaded regions 
correspond to time periods where wall effects prevent from 
proper fluorescence measurements. 



two side channels. A microfluidic switch (Upchurch) al- 
lows to exchange the solute solutions in the side channels, 
leading to a time-dependent tuning of the gradient. Par- 
ticles under investigation (colloids or A-DNA) are located 
in the center channel. We use polystyren-carboxylate flu- 
orescent colloids (F8888 200nm, Molecular Probes) in a 
Tris buffer (ImM, pH9), and DNA (48 kbp, Fermentas, 
Germany) labelled with YOYO-1 (Invitrogen) , in Tris- 
EDTA buffer (ImM, pH7.6). As a solute, we have con- 
sidered here either fluorescein or bare salts -LiCl, NaCl 
and KC1-. In a stationary regime, a linear profile of the 
solute concentration is achieved, as illustrated in Fig. [!]■ 
(b), with fluorescein as a benchmark. 

Let us first explore the response of colloidal particles 
under stationary salt gradients. Starting from a con- 
figuration with all particles gathered on one side of the 
channel, the salt gradient is switched toward the oppo- 
site channel wall: colloids are accordingly observed to 
drift toward the higher solute concentration, see Fig. [2} 
(a). In order to quantify this motion, we plot in Fig.[2]-(b) 
the time dependent location of the colloidal population, 
defined as the maximum of the distribution. The ob- 
served drift is close to linear, as would be expected for a 
constant mobility \i = Vdp/Vc. However a slight devia- 
tion from this linear expectation can be observed at long 
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22 ± 2 
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TABLE I: Values for the diffusiophoretic mobility Ddp 
extracted from Fig. [5] (see text). Experimental trap strength 
K extracted from the trapping experiments in Fig. [3] and 
compared to the predicted value from Eq. Q. 



times. This result can be understood by going more into 
details of the diffusiophoretic transport. Indeed, in the 
case of electrolytes as a driving solute, the mobility \i is 
expected to depend on the salt concentration c, scaling 
as /x(c) ~ c _1 , and the diffusio-phoretic velocity under a 
salt concentration gradient Vc(x, t) can be written as [3]: 



V DP = D D p V log c 



(1) 



with Ddp a diffusio-phoretic mobility pQ . Physically, this 
dependence originates in the balance between osmotic 
forces and visous stresses occuring within the Debye layer 
at the particle surface [2j [9] . This leads to Ddp scaling 
typically as Ddp ~ ^bT/t]£bi with r] the water viscosity 
and Ib the Bjerrum length (JLb — 0.7nm in water). 

This non-linear dependence on solute concentration is 
at the origin of the deviation from constant drift ob- 
served experimentally. Indeed under a linear salt profile 
c(x) = cq/2(1 ± 2x/£), with t~ x the slope of the concen- 
tration gradient (t~ x = Vc/co), and x the distance to 
mid-channel, Eq. ([!]) predicts that the position X of a 
colloid obeys dX /dt = ±D BP /(£/2 ± X ). This equa- 
tion can be solved analytically and provides a very good 
fitting expression for experimental data, see Fig. [2|(b), 
allowing for the determination of the colloids diffusio- 
phoretic mobility Dbp- In the case of colloids, the dif- 
fusiophoretic mobility was measured for three different 
salts, exhibiting salt specificity effects with values for the 
mobility in the order LiCl>NaCl>KCl, see Table I, in 
line with previous results [2]. Furthermore, we performed 
additional experiments with A-DNA in place of spherical 
colloids: Figs. [2]-(c,d) demonstrate a comparable motion 
of DNA molecules resulting from the diffusiophoretic mi- 
gration towards higher concentrations of salt, Table I. 

We now turn to the effect of solute time-dependent os- 
cillations: while in the previous section migration was 
studied under a stationary gradient, we now generate 
concentration gradients oscillations - Fig. [l]-(b) - (with 
period T in the range T /2 ~ 180 — 300s in the follow- 
ing). Typically, experiments show that a linear solute 
gradient establishes in the system within a time T t ~ 50s 
for £ = 800/im, in agreement with numerical estimates 
[T4] . The result is shown in Fig. ^ (a) for colloids: start- 
ing from a homogeneous distribution, the particles pop- 



3 



Colloids 




0.8 



0.6 



0.4 



0.2 




-g.i_ 0.1 ( 



FIG. 3: Trapping under oscillatory gradients: (Top): trap- 
ping of colloids, (a), and A-DNA, (b), under salt gradient 
oscillations. Buffers and salt concentration are identical to 
Fig. [2] with T = 600s (scale bar 100/xm). (Bottom): exper- 
imental particles density profile in stationary state, as ob- 
tained from the fluorescence images (symbols), fitted to a 
gaussian curve (dashed line). 



ulation evolves after a few cycles towards a band of gaus- 
sian shape with stationary width, Fig. [3]- (a) [14]. Note 
that the position of the band keeps oscillating around 
the center with the salt cycles, but its width remains sta- 
tionary. Moreover, the process is robust: starting from 
inhomogeneous particles distribution, or changing the os- 
cillation frequency yields the same width for the trapped 
band. Finally as for the colloids motility, the width of 
the trapped band is dependent on the salt nature: LiCl 
> NaCl > KC1 in terms of trapping efficiency. Addi- 
tionally, we also performed similar experiments with A- 
DNA as motile particles: Fig. |3]-(b) shows that trapping 
is equally achieved with macromolecules, the latter gath- 
ering also in a narrow band with gaussian profile. 

Going further, a physical interpretation for the segre- 
gation can be proposed: as we now show, trapping of 
the particles indeed results from the rectification of their 
motion under the oscillating driving field of the salt con- 
centration. The particles' population obeys the Smolu- 
chowski equation 

d t p = -V • (-D c V Po + D D pV[logc] x p ) , (2) 

coupled to the solute diffusive dynamics on c(x,t). Ne- 
glecting transients in the salt dynamics, its concentra- 
tion profile can be written as c(x,t) = co(l/2 + /(t) x/£), 
with f(t) a time-dependent function, oscillating with the 
forcing periodicity, To. Typically, f(t) can be approxi- 
mated as a ±1 step-function, so that (f(t)) = while 



(/(^) 2 ) = 1 7^ 0? w ^ n (') t ne time average over T . Ac- 
cordingly the position X (t) of the particle population 
is expected to oscillate with the gradients following the 
description above and Eq. ([!]). In the limit of small 
excursions around the center (Xq <C £), this reduces to 
dX /dt ~ 2D DP f(t)/£ [1 - 2f(t)(X /£) + ...]. 

Now, in the limit of fast salt oscillations, one expects 
the distribution in steady state to behave to leading order 
as p(x,t) ~ po(x — Xo(t)). Inserting this guess into the 
Schmoluchowski equation, Eq. ([2| and using the above 
equation for Xo(t) leads to the following equation for the 
particle density profile po : 



AD 



(J) = -D c \/ X p + ^^(/ 2 ) x Sx x po * (3) 

with Sx = x — Xo(t). In deriving this equation, we have 
furthermore averaged out over the fast salt variables. 
Solving this equation in the stationary state predicts a 
gaussian distribution for the particles: 
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that oscillates as a whole around the channel central posi- 
tion. The model thus reproduces the experimental obser- 
vations: trapping towards a gaussian distribution, with 
a frequency independent width. We furthermore assessed 
the validity of this description by performing a full nu- 
merical resolution of the coupled particle and solute dy- 
namics [T4] . 

Physically the origin of the focusing lies in the non- 
linear dependency of the diffusio-phoretic phenomenon 
versus the solute concentration, Eq. 0. Under a con- 
stant solute gradient, the velocity of particles is larger 
in regions with smaller solute concentration. Accord- 
ingly, the front particles move slower than the back ones: 
iterated over the oscillations, such a process leads to 
the observed focusing. The balance with Brownian mo- 
tion leads to an harmonic "osmotic" trapping potential, 
Vtrap(^) = \Kx 2 , with K = k^T/cr 2 . The experimental 
values for the measured trapping strength K are gath- 
ered in Table I for various salts. These values are in 
good agreement with theoretical predictions, in which 
T>dp was set to the mobility measured independently in 
channel-crossing experiments (Fig. [5| . A final important 
remark is that the trapping potential is indeed strong: 
over the scale of the system £, the trapping free-energy 
well has a depth of A J 7 « Kt 2 = k B T x ^gf up to 
hundreds of k B T ! (and independent of £). 



We now generalize the previous results to more complex 
geometries in order to demonstrate the robustness and 
generic character of the above scenario. To this end, we 
tested the localization phenomenon in a circular, "cell 
shaped", chamber, Fig. [4]- (a) -(b). Two different drivings 
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FIG. 4: Trapping of colloids in a circular chamber: (a)-(b) 
Initial distribution of the fluorescent colloids and sketch of 
the experimental set-up. The central circular well is 650/xm 
in diameter with gel walls 125 /xm wide (£ = 900/xm; scale bar 
200/im). The salt (LiCl) concentration oscillates in the two 
side channels (period To = 480s) either anti-symmetrically, 
Ci(t) - (Ci) = -[C 2 (t) - (C 2 )], or symmetrically d(t) = 
C2(t), between Co = lOOraM and (in addition to a TRIS 
buffer), (c)-(d): Stationary colloidal distribution under an 
antisymmetric driving (c) and under symmetric driving (d). 
(c , )-(d , )\ Theoretical predictions under corresponding exper- 
imental conditions. The predicted profiles are those obtained 
by avering out the oscillating salt distribution (see text). The 
values for Ddp were taken from experiments in Fig. [2] see 
Table I, while t and To take their experimental values. 

are explored: an antisymmetric driving, where the salt 
boundary concentrations, C\ and C2 are switched period- 
ically between and Co with antisymmetric phase; and a 
symmetric driving, where C\ = C2 switches periodically 
between and cq. As demonstrated in Figs. |4]-(c),(d), 
the solute oscillations again produce a localization of the 
particles in the cell chamber [13] . Furthermore the sym- 
metry of the pattern depends directly on the symmetry of 
the driving: linear - "cat eye" shape- for the antisymmet- 
ric driving, and circular -cell-shape- for the symmetric 
driving, thus demonstrating the versatility of the trap- 
ping process. A further remark is that this segregation 
process is also robust w.r.t. the initial distribution of the 
particles. 

A rationalization of the different patterns observed 
can be obtained along similar lines as above. First 
the solute diffusion equation is solved with the ap- 
propriate boundary conditions under periodic inver- 
sion. In the case of the antisymmetric geometry, the 
salt concentration profile c(r, 6,t) (with r, the ortho- 
radial coordinates) is found to take the general form 
c(r,0,t) = cq/2[1 + 0)f(t)] With (JantiM) = 

n odd ( r /^) /C sin[fc0], and f(t) the periodic step- 
like function with period Tq. In the case of a symmet- 



ric driving, one finds c(r, Q,t) = cq/2 [1 + £ sym (r, 0, t)] 
where 5 8ym (r,t) = ± £ fc>odd (l/fc)Im [/fc( r ) ex PQ' k uj *)] 
with f k (r) = I (r/5 k )/I (e/25 k ), 5 k = ^jD s /kuj the 
salt diffusive length, uo = 2tt/To and Io(x) is the Bessel 
function of order 0. In both cases, averaging out the 
fast salt variables, the trapping potential is then ob- 
tained in the form Vt raiP (r,0) = ^(-Ddp/^c) x <K r ?#) : 
0(r, 0) = £ ant i(r, #) 2 for the antisymmetric case, while 
(j){r,0) ~ Kor 4 for the symmetric driving, with 1 = 
Wndl cosh(l/\/25i)/l. These predicted localization pat- 
terns are exhibited in Figs. Q-(c , ),(d / ), showing a good 
qualitative and semi-quantitative agreement with the ex- 
perimental results. 

To conclude, we have demonstrated a new mecha- 
nism leading to segregation and pattern formation of col- 
loids and macromolecules, originating in rectified time- 
dependent solute contrasts. A key ingredient of the pro- 
cess is the so-called diffusiophoresis, a passive transport 
phenomenon leading to motion of particles under chem- 
ical potential gradients. Our results highlight the im- 
portance of solute contrast induced transport in out-of- 
equibrium processes, with potential implications in soft 
matter, as well as in living and chemical systems, where 
concentration gradients are ubiquituous. 
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